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There has been a growing interest in constructing exact ground states for correlated 
Fermion systems in certain parameter regions [1-7]. For one dimensional quantum systems, 
an effective method to construct ground states is the so called matrix product approach JT| . 
Its origin can be traced back to the S — 1 spin chain model 0. For higher dimensional 
correlated Fermion systems, there are also some approaches for the construction of exact 
ground states ||. In j| the conditions for 77-pairing, paramagnetic, Neel, charge-density- 
wave and ferromagnetic states to be exact ground states of the generalized Hubbard model 
have been presented by using the so called optimal ground state approach. These approaches 
recovered and improved the results obtained previously. 

In this paper we give some new exact ground states for the generalized Hubbard models. 
We first describe a general way to construct some exact ground states for homogeneous 
quantum systems. Let H be the Hamiltonian of a quantum lattice system (correlated 
Fermion system) of the form H = J2<ij> hij, where < ij > denotes the nearest neighbours, 
i, j are points on a finite lattice T with A lattice sites, H acts in the tensor space 7i associated 
with T (i.e. 7i is the A-fold power of some given Hilbert space K associated with a single 
lattice point), hij is a self-adjoint operator acting on the tensor product K <S> K associated 
with We call the system homogeneous if the spectral decomposition of is independent 
of For square lattice case, this homogeneity implies that the horizontal and vertical 
couplings are the same. The exact ground states of H are by definition the eigenstates 
belonging to the lowest eigenvalue of H. Let |f Q >, a = 1, ...,n be n possible states of a 
local site, for instance, empty, spin up (down) and double occupied states. We have, for 
some given indices a, /3 6 1, ...,n, the following conclusions: 

I) If |£ Q >i \£ a >j is an eigenvector of and the corresponding eigenvalue is the lowest 
one, the following state is a ground state of the Hamiltonian H of the system, 

^ = n i*« >< • w 

i=i 

II) If \C a >i |f a >j, |f/3 >i If/3 >j and |£ Q >i 1^ >j +|f/3 >i \£ a >j are eigenvectors of 
with the same and lowest energy, the state 

^ U = I[(l£a >i +lf/3 >i) (2) 

1=1 

is a ground state of H. 
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III) If |£ Q >i \£ a >j and |£ Q >; \£p >j >j |£ a >^ are eigenvectors of ^ with the 
same and lowest energy, the state 

tf" = J2(-l) i+1 \Ca >1 |6* >2 ..-1^ >i -16* >A (3) 
i=l 

is a ground state of if. 

IV) If |£ a >j |£ Q >j and |£ a >j >j |£ Q >j are eigenvectors of /iy with the 
same and lowest energy, the state 

i> IV = E 16* >1 16* >2 >i -|6* >A (4) 

i=l 

is a ground state of H. 

[Proof]. I is a consequence of homogeneity, II holds by a simple calculation. As for III 
(IV can be similarly proved) we proceed as follows. Let 

16* > l o 

16* >i +163 >i H6* >i 

We have 

/ 16, >i 16* >j 
A- A = I 

V |6* >i |f/J >j -16? >^ 16* >j 16* >* 16* >i 

The nonzero elements of AjA,- satisfy h^a^ = eaij, with e the lowest eigenvalue of hij. 
Therefore the elements of the matrix M = n^=i A, as well as their linear combinations are 
ground states of H with the same lowest eigenvalue Pe, where P is the number of nearest- 
neighbour pairs of the lattice. The linear combination of the elements M(i^ and M( 2 ,i) gives 
rise to the ground state ■ 

The state ip n is a sum of terms containing factors of the form |£ Q >j and |^ >_,-. If 
there exists a local operator o such that o|£ Q >= ai|£ Q >, >= 02)63 >, ai 7^ 02 and the 
operator O = X)f 0j commutes with the Hamiltonian if, then 

A 

^ = E 16* >1 -|& >mi -I6? >™ ; -16, >A (5) 
mi,...,m ( =l 

are also ground states for all / = 0, 1, A. In this case the states V^i^A-i (resp. V'fio.A) are 
just the states V ;/y (resp. ■?//). However, the state ?/> /y (resp. has a larger parameter 
region than f// 7 . In addition, all the states xjjf 1 are degenerate. 
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We now consider the generalized Hubbard model on an arbitrary d- dimensional lattice 
with Hamiltonian ||, 

+1 (fat - I) (mi - I) + ("it - I) - I)) + ^ - " !) 

(6) 

where cf CT and c^, cr =|, j, are canonical Fermi creation resp. annihilation operators, nj )(7 = 
cl^Cjo- and rij = + riij, are the number operators, = cj+csj, S 1 " = c|, c^-p , 5? = (n»f — 
rij|)/2; is the chemical potential and Z the coordination number of the <i-dimensional 
lattice; t, X, U, V, Y, J xy and J z are real valued coupling constants for, respectively, single 
particle hopping, bond-charge interaction, on-site Coulomb interaction, nearest-neighbour 
Coulomb interaction, pair-hopping and XXZ-type spin interactions. 

The well known states related to the Hubbard model are the //-pairing, paramagnetic, 
Neel, charge-density-wave (CDW) and ferromagnetic (F) states, defined respectively as fol- 
lows: 

\Neel >= UiEB 4 T IW 4J0 >, >= (eJ =1 e ir »c] L c}T |0 >, m e N 

\para >= l\ ieA cJ T ft ieA , cj ; |0 >, |CLW >= ILes cjfcJjO >, |F >= Hi 4^0 >, 

where |0 > is the vacuum state, A n A' = 0, A and A' together span the whole lattice, B 
and B 1 are the odd and even sublattices of a bipartite lattice, rj is a real parameter. The 
^-pairing states display off diagonal long-range order (ODLRO) M, i.e. superconductivity 
(Meissner effect and flux quantization |J). 

The exact ground states of H depend on the coupling constants. At certain parameter 
regions the above states (0) become exact ground states 0, |4[], as also follows from our 
conclusions. We observe that some of the ground states could be physically equivalent 
owing to some symmetries of the system in certain parameter regions, for instance, when 
Jz = Jxy, H is SU(2) invariant. Nevertheless besides the well known exact ground states 
the quantum system ([]) has many additional exact ground states of types ip n , if} 111 and ip IV 
with respect to different parameter regions. In the following we simply list some of these 
new ground states. 



4 



We first present some exact ground eigenstates of the form ifj 111 . 

M >= ^E(-l) m cLcL-4_ CT ...cL|0 >, N = (8) 

1=1 

is a ground state when t = X, J z = -J xy < 0, -| > max(^f + ^r — V, 4f + ^ + 2\t\, + 
V + Y, f - f + 2|t|, f - f - V, ± + V-Y). The eigenvalue is (-§ + f + f )P. If 
one sets J 2 = J xy = instead of J 2 = — J xy < 0, state @ is still a ground state but is then 
equivalent to the \para > state, because > is a linear combination of \para > states, 
while under condition t — X , J z — J xy = 0, any \para > state is an exact ground state. 

The state 

M >= E(-i) i+1 i4-4r c l-cL|o > (9) 

1=1 

is a ground state at parameter region X + ^<0ift>0, X + §<tift<0, J z = 2 jf + ^± - 
4(t-2X) andmax(-t+2X + f -V, -t+2X+V+Y+%, -t+2X-V-%) <^< V -Y + 
f -t+2X, If > max(^-t+2X+§ , t-2X-± -£+2(X-t) 2 /(V-Y-t+2X+% -£)), or 
g > max(-t+2X+f -V, -t+2X+\/+F+f , -t+2X-V-%), £ > V-Y + ±-t+2X, 
t-2X-%-^+2(X-t) 2 /(V-Y-t+2X+%-§) > J -f> ^-t+2X+§ for t > 2X, t ^ X, 
or max(t-V+f, t+^+V+Y, t-V-%, t+V+^-Y) <£< min{ J -f-t-^, -^f-t-f), 
J z = 2U + M + 4t ? ^ < o for t = X < 0. The eigenvalue of |V>f > is (2X - t + f )P. 

For £ < and J 2 = 4(t + ^ - |) the following state 

HI >= Nj2(-iy+ 1 c{A,„---h...c{,Jo >, (io) 

i=\ 

1 being the identity operator, is an exact ground state, for t ^ X, -| > X > t — & and 
V-Y + t-% >§> max(t+^-V, t-^ + V + Y, t-V-f), J -f- >max(£+t-%, f- 
t-£ + 2(X-t) 2 /(V-Y + t-£-$)) or ^ > max{t + §-V, t- § + V + Y, t-V- f ), 

JL>y-y + *-!, |-t-^ + 2 (x-t) 2 /(y-y + t-^-|)>^>§ + i-|, 

or for t = X < 0, /i > and max{t + | - V, t-^ + V + Y, t-^ + V-Y)<^< 

Another kind of exact ground states of the generalized Hubbard model of the form 
is given by: 

m >= E c UL-4, T 4,i- c Llo > . (ii) 
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The corresponding parameter region is, fort < 2X, t ^ X , J z = 4(|g+#+t— 2X), X— ^ > t 
if t > 0, X-§ > if t < 0, ^ > max{t-2X-V+^ , t-2X + V + Y+%, t-2X-%-V) 
and g < V-F + t-2X + |, > max(f + 2t-4X + f , 4X-2t- f - f + 4(X-t) 2 /(V- 

Y + t-2X + or ^ > ^-F + t-2X + f, 4X - 2t - § - f + 4(X - 1) 2 /(V - Y + 
t - 2X - g + f ) > > f + 2t - 4X + f , and for t = X > 0, J z = 4(g + § - t), Li < 0, 
mm(% + f-f,t-^-f) > § > max(-t + V + Y + §, -t - f -V, -t + V-Y + %). 
\tp% > corresponds to the eigenvalue (t — 2X + %)P. 

For t > and J 2 = 4(-t + ^ - §) the state 

A 

|^ >= X^c^ct-li-c^lO > (12) 
t=i 

is an exact ground state, for t X, —% < X < and V— F — £— £ > ^ > maa;(-t+|- 
V, -t-$+V+Y, -t-V-f), *f. > m ax{^-t-^, ± +t -± + 2{X-tf / {V -Y -t-±-± )) 
or § > max{-t + § - V, -t-f + V + y, -t-V-f), JL>V-y_t-|, 
f + 1 - ^ + 2(X - t) 2 /(V - Y - 1 - £ - § ) > % > § - 1 - § , or for t = X, // > 0, and 
mox(-t+|-V; -t-f + V + F, -t-^ + V-Y) < £ < rrtin{t+% + !f , t-?f + %. 
The eigenvalue is (— t + #)P. 

The states of the form f/> 7/ are sums of states with different number of electrons. As the 
Hamiltonian commutes with the number operator, any projection of a ground state of the 
type ip 11 to a state with a given number of electrons is also a ground state. These ground 
states are of the form (|5|), for instance, for I = 2, A — 2, the following states are ground 
states, 

A 

|V>f>= E c\ r c{ l ..J mia ...cl ia ..J M c{ l \0> (13) 

mi,...,mi=l 

if the parameters satisfy J z = 4(^ - t + §), V = -§ - t - ^, y. < 0, \Y\ - 2t < 
| < 2t - %^ - f for £ = X > 0. For t ^ X we have J 2 = 4(§ + t - 2X + § ), 

V = -J±+t-2X-£, Li < 0, X > max(f ,£ + §,§) and 2t - 4X + Y < § < 2t - 4X - y, 

> max(f +2£-4X+f , 4X-2t-f -§+4(X-t) 2 /(2t-4X-y-§)) or 2£-4X+y < f , 
2i-4X-y < f , 4X-2t-f -| + 4(X-t) 2 /(2t-4X-y-f) > 3 -f- > | + 2t-4X + f . 
The eigenvalues are (^ + t — 2X)P. For / = 0, 1, A — 1, A, the states > are already 
included in ip 1 and ip IV . 
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Similarly we have the following ground states: 



A 



> £ 4 lCT ...4 iCT |0>, / = 2,..,A-2 (14) 

mi,...,rri(=l 



when 



and, for t ^ X, t > 0, K - 2t < § < -F - 2t, n > 0, § + t > X > -f , > 
max(-t + - § + t+ § - 2(X - t) 2 /(§ + Y + 2t)) or F - 2t < § , -Y - 2t < §, 
M >0, f +t>X>-|, + -2(X-t) 2 /(f + F + 2t)>^>-t-f + ^,or 
for t = X > 0, 

-2t+|r|<|<2t-^ + ^, „>0. (16) 

The corresponding eigenvalue is — t)P. 

From formulae (1-4) one can construct more exact ground states such as >= 
^Ef =1 (-l) i+1 c t ljt c t u 4 jT 4 ii ... C t J ...4 >T cJ Vji iO >, \r 7 >= ArEf =1 (-l) i+1 cJ i(T |0 >, \tf >= 

NJ2f=i c i,t c U c 2,T c 2,|"- c iff-" c A,T c A,|l^ > ' 1^9 >= ^SiLi4,<r|0 > i n certain parameter re- 
gions. Some of the states have the same eigenvalues. For instance, the state > has 
the same eigenvalue as >, however |-0f > and > correspond to different parameter 
regions. Thus in general the eigenvalues are not degenerate. 

The hij in (^) have two different sets of eigenvectors and eigenvalues for t — X and 
t ^ X respectively. Therefore the parameter regions of the ground states have a jump 
when t — > X. In addition, all the ground states we presented above have no ODLRO. 
A ground state with ODLRO can be obtained from formula ([D by taking |£ a >= |0 > 
and >= cjcj|0 >. One then gets ground states that are equivalent to the ^-pairing 

states with 77 = 0: |^ =0 >= Emi,...,™^ c L 1 T c mi|-- c m ! T c m ! |l° > within the parameter region 
§ < min{-2\t\ - 2V, f - V, - V - ^ - f), t = X, Y = 2V, /1 = and V < 0. 

Some properties of the phase diagram of the ground states can be analysed according to 
their related parameter regions. For instance, taking into account the conditions ([15]) and 
flTfiD (for simplicity we set J xy = 0), we have that iV'f -0 > is a ground state for ^ < \2t\. 
When § < -2t, the 77-pairing > is no longer a ground state, but |02« > given by (JT3J) 

is a ground state. 
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Some relations between the ground states above and the ground states studied in j|] can 
be discussed according to their parameter regions. As an example we set t = X, Y = fi = 0, 
V > 0, and J z = —J xy < 0. It can be shown that for \t\ > J xy /8 the states > and \F > 
are ground states in regions I and IV. While the state \CDW > is a ground state in region II 
and III, see Fig. 1. More properties can be obtained when one changes the relations among 
the parameters. For instance, if t = V, then the states > and \F > (resp. \CDW >) 
are ground states only in region I (resp. region II). 
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Figure 1: A ground state phase diagram related to states \F >, \CDW > and \ip° > 
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